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QCD in a strong magnetic field yields an example of a rich, sophisticated and 
controllable dynamics. 



1. Prologue 

On the first day of this wonderful meeting, I decided to change the topic 
of my talk. Already during the first Session, I felt it necessary to present 
something special, something connected with Arkady. So, I decided to 
talk about the dynamics of QCD in a strong magnetic field. This talk, it 
seemed to me, was appropriate indeed. The point is that the last time when 
I interacted with Arkady in person was in April 1994, in the Institute for 
Theoretical Physics in Santa Barbara. At that time, we, Valery Gusynin, 
Igor Shovkovy, and myself, had just finished our first work [ 1] in the series 
of papers concerning the role of a magnetic field in dynamical symmetry 
breaking in 2+1 and 3+1 dimensional field theories (actually, at that time, 
we discussed only the 2+1 dimensional case). So, it seemed to be appropri- 
ate to "continue" that discussion at the Arkadyfest in Minneapolis 8 years 
later. 

I was lucky with this decision. At the time of its presentation, this work 
had not been yet completed. During my talk, Arkady raised a question 
concerning the dynamics of the magnetic catalysis in QCD with a large 
number of colors, N c — > oo. The answer to this question, Igor Shovkovy 
and I got the next day, helped to finish the work [ 2] . 



2. Introduction 

Since the dynamics of QCD is extremely rich and complicated, it is impor- 
tant to study this theory under external conditions which provide a control- 
lable dynamics. On the one hand, this allows one to understand better the 
vacuum structure and Green's functions of QCD, and, on the other hand, 
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there can exist interesting applications of such models in themselves. The 
well known examples are hot QCD and QCD with a large baryon density. 

Studies of QCD in external electromagnetic fields had started long ago 
[3,4]. A particularly interesting case is an external magnetic field. Us- 
ing the Nambu-Jona-Lasinio (NJL) model as a low energy effective theory 
for QCD, it was shown that a magnetic field enhances the spontaneous 
chiral symmetry breakdown. The understanding of this phenomenon had 
remained obscure until a universal role of a magnetic field as a catalyst of 
chiral symmetry breaking was established in Refs. [1,5]. The general result 
states that a constant magnetic field leads to the generation of a fermion 
dynamical mass (i.e., a gap in the one-particle energy spectrum) even at 
the weakest attractive interaction between fermions. For this reason, this 
phenomenon was called the magnetic catalysis. The essence of the effect is 
the dimensional reduction D — > D — 2 in the dynamics of fermion pairing 
in a magnetic field. In the particular case of weak coupling, this dynamics 
is dominated by the lowest Landau level (LLL) which is essentially D — 2 
dimensional [1,5]. 

The phenomenon of the magnetic catalysis was studied in gauge theo- 
ries, in particular, in QED [ 6,7] and in QCD [ 8,9]. In the recent work [ 9], 
it has been suggested that the dynamics underlying the magnetic catalysis 
in QCD is weakly coupled at sufficiently large magnetic fields. Here we will 
consider this dynamical problem rigorously, from first principles. In fact, 
we show that, at sufficiently strong magnetic fields, \eB\ » Ag C£) , there 
exists a consistent truncation of the Schwinger-Dyson (gap) equation which 
leads to a reliable asymptotic expression for the quark mass m q . Its explicit 
form reads: 



m, 



~2Ci|e gJ B|(c g a s ) 2/3 exp 



4N r TT 



a s (N2-l)\n(C 2 /c q a s ) 



(1) 



where e q is the electric charge of the q-th quark and N c is the number of col- 
ors. The numerical factors C\ and Ci equal 1 in the leading approximation 
that we use. Their value, however, can change beyond this approximation 
and we can only say that they are of order 1. The constant c q is defined as 
follows: 



(2N U + N d ) 

D7T 



(2) 



where N u and Nd are the numbers of up and down quark flavors, respec- 
tively. The total number of quark flavors is Nf = N u + Nd- The strong 
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coupling a s in the last equation is related to the scale ^/|e_B|, i.e., 



1 , , \eB\ , , UN C -2N f 

— ~61n-^— -, where b= — '-. (3) 

A* CD ' 127T 1 ' 



We should note that in the leading approximation the energy scale \/\eB\ 
in Eq. (3) is fixed only up to a factor of order 1. 

The central dynamical issue underlying this dynamics is the effect of 
screening of the gluon interactions in a magnetic field in the region of mo- 
menta relevant for the chiral symmetry breaking dynamics, nil <C \k 2 \ <C 
\eB\. In this region, gluons acquire a mass M g of order y/N f a s \e q B\. This 
allows to separate the dynamics of the magnetic catalysis from that of con- 
finement. 

Since the background magnetic field breaks explicitly the global chiral 
symmetry that interchanges the up and down quark flavors, the chiral sym- 
metry in this problem is SU(N U ) L x SU(N U ) R x SU(N d ) L x SU(N d ) R x 
U(~\1)a- The U(~\V)a is connected with the current which is an anomaly 
free linear combination of the U^(1)a an d U^ u \1)a currents. The gen- 
eration of quark masses breaks this symmetry spontaneously down to 
SU(N U ) V x SU(N d ) v and, as a result, + TVj - 1 gapless Nambu- 
Goldstone (NG) bosons occur. In Sec. 4, we derive the effective action 
for the NG bosons and calculate their decay constants and velocities. 

The present analysis is heavily based on the analysis of the magnetic 
catalysis in QED done by Gusynin, Miransky, and Shovkovy [ 7] . A crucial 
difference is of course the property of asymptotic freedom and confinement 
in QCD. In connection with that, our second major result is the derivation 
of the low energy effective action for gluons in QCD in a strong magnetic 
field. The characteristic feature of this action is its anisotropic dynamics. In 
particular, the strength of static (Coulomb like) forces along the direction 
parallel to the magnetic field is much larger than that in the transverse 
directions. Also, the confinement scale in this theory is much less than that 
in QCD without a magnetic field. This features imply a rich and unusual 
spectrum of light glueballs in this theory. 

A special and interesting case is QCD with a large number of colors, in 
particular, with N c — > oo (the 't Hooft limit). The question about it was 
raised by Arkady during my talk. This theory is considered in Sec. 6. 
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3. Magnetic catalysis in QCD 

We begin by considering the Schwinger-Dyson (gap) equation for the quark 
propagator. It has the following form: 



G~ 1 (x,y) = S- 1 (x,y) + 4na s ^ / G(x, z)T v (z, y, «)£>„„(«, x)d i zd i u{A) 



where S(x, y) and G(x, y) are the bare and full fermion propagators in an ex- 
ternal magnetic field, D vtl (x, y) is the full gluon propagator and T v (x, y, z) 
is the full amputated vertex function. Since the coupling a s related to the 
scale \eB\ is small, one might think that the rainbow (ladder) approxima- 
tion is reliable in this problem. However, this is not the case. Because of 
the (l+l)-dimcnsional form of the fermion propagator in the LLL approx- 
imation, there are relevant higher order contributions [6,7]. Fortunately 
one can solve this problem. First of all, an important feature of the quark- 
antiquark pairing dynamics in QCD in a strong magnetic field is that this 
dynamics is essentially abelian. This feature is provided by the form of 
the polarization operator of gluons in this theory. The point is that the 
dynamics of the quark-antiquark pairing is mainly induced in the region 
of momenta k much less than ^/|e_S|. This implies that the magnetic field 
yields a dynamical ultraviolet cutoff in this problem. On the other hand, 
while the contribution of (electrically neutral) gluons and ghosts in the 
polarization operator is proportional to k 2 , the fermion contribution is pro- 
portional to \e q B\ [ 7]. As a result, the fermion contribution dominates in 
the relevant region with k 2 <C \eB\. 

This observation implies that there are three, dynamically very different, 
scale regions in this problem. The first one is the region with the energy 
scale above the magnetic scale ^/|e£>|. In that region, the dynamics is 
essentially the same as in QCD without a magnetic field. In particular, 
the running coupling decreases logarithmically with increasing the energy 
scale there. The second region is that with the energy scale below the 
magnetic scale but much larger than the dynamical mass m q . In this region, 
the dynamics is abelian like and, therefore, the dynamics of the magnetic 
catalysis is similar to that in QED in a magnetic field. At last, the third 
region is the region with the energy scale less than the gap. In this region, 
quarks decouple and a confinement dynamics for gluons is realized. 

Let us first consider the intermediate region relevant for the magnetic 
catalysis. As was indicated above, the important ingredient of this dy- 
namics is a large contribution of fermions to the polarization operator. It 
is large because of an (essentially) 1+1 dimensional form of the fermion 
propagator in a strong magnetic field. Its explicit form can be obtained 
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by modifying appropriately the expression for the polarization operator in 
QED in a magnetic field [ 7] : 

N f 

q—l Q 

N f 

pAB,jxv ^ _C± sAB ^ _ ^ £ 1^1 , m » <<; |fcf| <<; |ejB|; (6) 

5=1 I 

where <?||" y = diag(l, 0, 0, — 1) is the projector onto the longitudinal sub- 
space, and fcy* = g^"k v (the magnetic field is in the x 3 direction). Similarly, 
we introduce the orthogonal projector g^ 1 = g^ v — g^ u = diag(0, —1,-1,0) 
and = g^"k v that we shall use below. Notice that quarks in a strong 
magnetic field do not couple to the transverse subspace spanned by g^ 
and fc^ . This is because in a strong magnetic field only the quark from the 
LLL matter and they couple only to the longitudinal components of the 
gluon field. The latter property follows from the fact that spins of the LLL 
quarks are polarized along the magnetic field (see the second paper in [5]). 

The expressions (5) and (6) coincide with those for the polarization 
operator in the 1 + f dimensional massive QED (massive Schwinger model) 
[ 10] if the parameter a s \e q B\/2 here is replaced by the dimensional coupling 
ot\ of QED\ + \. As in the Schwinger model, Eq. (6) implies that there is a 
massive resonance in the k^kY, — k?,gf^ v component of the gluon propagator. 
Its mass is 

N f 

M* = J2^\e q B\ = (2N u + N d )^-\eB\. (7) 

9=1 

This is reminiscent of the pseudo-Higgs effect in the (l+I)-dimensional 
massive QED. It is not the genuine Higgs effect because there is no complete 
screening of the color charge in the infrared region with <C m 2 q . This 
can be seen clearly from Eq. (5). Nevertheless, the pseudo-Higgs effect is 
manifested in creating a massive resonance and this resonance provides the 
dominant forces leading to chiral symmetry breaking. 

Now, after the abelian like structure of the dynamics in this problem 
is established, we can use the results of the analysis in QED in a mag- 
netic field [ 7] by introducing appropriate modifications. The main points 
of the analysis are: (i) the so called improved rainbow approximation is 
reliable in this problem provided a special non-local gauge is used in the 
analysis, and (ii) for a small coupling a s (a in QED), the relevant region 
of momenta in this problem is rr? q <§C |fc 2 | -C |e£?|. We recall that in the 
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improved rainbow approximation the vertex T v (x,y, z) is taken to be bare 
and the gluon propagator is taken in the one-loop approximation. More- 
over, as we argued above, in this intermediate region of momenta, only 
the contribution of quarks to the gluon polarization tensor (6) matters. [It 
is appropriate to call this approximation the "strong-field-loop improved 
rainbow approximation" . It is an analog of the hard-dense-loop improved 
rainbow approximation in QCD with a nonzero baryon density]. As to the 
modifications, they are purely kinematic: the overall coupling constant in 
the gap equation a and the dimensionless combination M 2 /\eB\ in QED 
have to be replaced by a s (N 2 — 1)/2N C and Mg/\e q B\, respectively. This 
leads us to the expression (1) for the dynamical gap. 

After expressing the magnetic field in terms of the running coupling, 
the result for the dynamical mass takes the following convenient form: 



m„ 



2d 



* 'z i \2/3 
A QCD \CqOis) exp 



1 

ba s 



a s (N 2 -l)\n(C 2 /c q a a ) 



■ (8) 



As is easy to check, the dynamical mass of the u-quark is considerably larger 
than that of the d-quark. It is also noticeable that the values of the u-quark 
dynamical mass becomes comparable to the vacuum value m^ n ~ 300 MeV 
only when the coupling constant gets as small as 0.05. 

Now, by trading the coupling constant for the magnetic field scale \eB\, 
we get the dependence of the dynamical mass on the value of the external 
field. The numerical results are presented in Fig. 1 [we used C\ — C2 — 
1 in Eq. (8)]. As one can sec in Fig. 1, the value of the quark gap 
in a wide window of strong magnetic fields, A 



QCD 



< \eB\ < (10 TeV) 2 , 
300 MeV in 



remains smaller than the dynamical mass of quarks "m dyn 
QCD without a magnetic field. In other words, the chiral condensate is 
partially suppressed for those values of a magnetic field. The explanation 
of this, rather unexpected, result is actually simple. The magnetic field 
leads to the mass M g (7) for gluons. In a strong enough magnetic field, 
this mass becomes larger than the characteristic gap A in QCD without 
a magnetic field (A, playing the role of a gluon mass, can be estimated 
as a few times larger than Kqcd)- This, along with the property of the 
asymptotic freedom (i.e., the fact that a s decreases with increasing the 
magnetic field), leads to the suppression of the chiral condensate. 

This point also explains why our result for the gap is so different from 
that in the NJL model in a magnetic field [3]. Recall that, in the NJL 
model, the gap logarithmically (i.e., much faster than in the present case) 
grows with a magnetic field. This is the related to the assumption that both 
the dimensional coupling constant G = g/A 2 (with A playing a role similar 
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Figure 1. The dynamical masses of quarks as functions of \n(\eB\/ Aq CD ) for N c = 3 
and two different values of Nf = N u + N d - (i) masses of u-quark (solid line) and d-quark 
(dash-dotted line) for N u = 1 and = 2; (ii) masses of M-quark (dashed line) and d- 
quark (dotted line) for N u = 2 and Nd = 2. The result may not be reliable in the weak 
magnetic field region (shaded) where some of the approximations break. The values of 
masses are given in units of Aqcd = 250 MeV. 



to that of the gluon mass in QCD), as well as the scale A do not dependent 
on the value of the magnetic field. Therefore, in that model, in a strong 
enough magnetic field, the value of the chiral condensate is overestimated. 

The picture which emerges from this discussion is the following. For 
values of a magnetic field \eB\ < A 2 the dynamics in QCD should be 
qualitatively similar to that in the NJL model. For strong values of the 
field, however, it is essentially different, as was described above. This in 
turn suggests that there should exist an intermediate region of fields where 
the dynamical masses of quarks decreases with increasing the background 
magnetic field. 

4. Effective action of NG bosons 

The presence of the background magnetic field breaks explicitly the global 
chiral symmetry that interchanges the up and down quark flavors. This 
is related to the fact that the electric charges of the two sets of quarks 
are different. However, the magnetic field does not break the global chi- 
ral symmetry of the action completely. In particular, in the model with 
the N u up quark flavors and the Nd down quark flavors, the action is 
invariant under the chiral symmetry SU(N U ) L x SU(N U ) R x SU(N d ) L x 
SU(N d ) R x £/(~)(1)a. The U^\1)a is connected with the current which is 
an anomaly free linear combination of the U^(1)a and U^ u \1)a currents. 
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[The U(~\1)a symmetry is of course absent if either N d or N u is equal to 
zero] . 

The global chiral symmetry of the action is broken spontaneously down 
to the diagonal subgroup SU(N u )v x SU(N d )v when dynamical masses of 
quarks are generated. In agreement with the Goldstone theorem, this leads 
to the appearance of AT^+iVj — 1 number of the NG gapless excitations in the 
low-energy spectrum of QCD in a strong magnetic field. Notice that there 
is also a pseudo-NG boson connected with the conventional (anomalous) 
U(1)a symmetry which can be rather light in a sufficiently strong magnetic 
field. 

Now, in the chiral limit, the general structure of the low energy action for 
the NG bosons could be easily established from the symmetry arguments 
alone. First of all, such an action should be invariant with respect to 
the space-time symmetry 5*0(1, 1) x SO (2) which is left unbroken by the 
background magnetic field [here the S0(1,1) and the S0(2) are connected 
with Lorentz boosts in the xq — X3 hyperplane and rotations in the x\ — X2 
plane, respectively]. Besides that, the low-energy action should respect the 
original chiral symmetry SU(N U ) L x SU(N U ) R x SU(N d ) L x SU(N d ) R x 
U^~\1)a- These requirements lead to the following general form of the 
action: 

L NG ~ 4 tr (g^d^ud^l + vlg^d^d^l 
+ f -ftv (ffp^&Et + v 2 d g^d^ d d^ d ) 

+ ^ (fffdntdjf + v 2 g^d^d^) . (9) 

The unitary matrix fields S u = exp (i J2a=i * ^ Alr u / fu) , = 
exp (iY^A^i 1 ^ Alr d /.fdj, and E = exp (i^jr/f^J describe the NG bosons 
in the up, down, and U^~\1)a sectors of the original theory. The decay 
constants /„, f d , f and transverse velocities v u ,v d ,v can be calculated by 
using the standard field theory formalism. The results for the + iVj — 2 
NG bosons in the up and down sectors, assigned to the adjoint representa- 
tion of the SU(N U ) V x SU(N d ) v symmetry, are [ 2] 

fu = §i\*B\, (10) 

The remarkable fact is that the decay constants are nonzero even in the 
limit when the dynamical masses of quarks approach zero. The reason of 
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that is the 1 + 1 dimensional character of this dynamics. A similar situation 
takes place in color superconductivity: in that case the 1 + 1 dimensional 
character of the dynamics is provided by the Fermi surface. 

Notice that the transverse velocities of the NG bosons are equal to 
zero. This is also a consequence of the 1 + 1 dimensional structure of the 
quark propagator in the LLL approximation. The point is that quarks 
can move in the transverse directions only by hopping to higher Landau 
levels. Taking into account higher Landau levels would lead to nonzero 
velocities suppressed by powers of |m g | 2 /|eS|. In fact, the explicit form of 
the velocities was derived in the weakly coupled NJL model in an external 
magnetic field [see Eq. (65) in the second paper of Ref. [ 5]]. It is 

2 l"Vd| 2 , \eB\ 

^~Tsr ln K^ <<L (12) 

A similar expression should take place also for the transverse velocities of 
the NG bosons in QCD. 

The decay constant / of the singlet NG boson connected with the spon- 
taneous breakdown of the U^~\1)a is [ 2] 

f = Nj = Y2^N] |eS| - (13) 

Its transverse velocity is of course zero in the LLL approximation. 



5. Anisotropic confinement of gluons 

Let us now turn to the infrared region with |fc| < md, where all quarks de- 
couple (notice that we take here the smaller mass of d quarks). In that re- 
gion, a pure gluodynamics realizes. However, its dynamics is quite unusual. 
The point is that although gluons are electrically neutral, their dynamics 
is strongly influenced by an external magnetic field, as one can see from 
expression (5) for their polarization operator. In a more formal language, 
while quarks decouple and do not contribute into the equations of the renor- 
malization group in that infrared region, their dynamics strongly influence 
the boundary (matching) conditions for those equations at k ~ m^. A con- 
ventional way to describe this dynamics is the method of the low energy 
effective action. This low effective action was derived in Ref. [ 2] . Here we 
will discuss its main properties. 

The low energy effective action is relevant for momenta |fc| < m^. Notice 
the following important feature of the action: its "bare" coupling constant 
g, related to the scale m<j, coincides with the value of the vacuum QCD 
coupling related to the scale \/\eB\ (and not to the scale ma)- This is 
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because g is determined from the matching condition at |fc| ~ m^, the 
lower border of the intermediate region rrid < \k\ < \/\eB\, where, because 
of the pseudo-Higgs effect, the running of the coupling is essentially frozen. 
Therefore the "bare" coupling g indeed coincides with the value of the 
vacuum QCD coupling related to the scale ^/|eB|: g = g s . Since this value 
is much less that that of the vacuum QCD coupling related to the scale mj, 
this implies that the confinement scale Xqc d of the action should be much 
less than Aqcd in QCD without a magnetic field. 

Actually, this consideration somewhat simplifies the real situation. 
Since the LLL quarks couple to the longitudinal components of the po- 
larization operator, only the effective coupling connected with longitudinal 
gluons is frozen. For transverse gluons, there should be a logarithmic run- 
ning of their effective coupling. It is clear, however, that this running 
should be quite different from that in the vacuum QCD. The point is that 
the time like gluons are now massive and their contribution in the running 
in the intermediate region is severely reduced. On the other hand, because 
of their negative norm, just the time like gluons are the major players in 
producing the antiscreening running in QCD (at least in covariant gauges). 
Since now they effectively decouple, the running of the effective coupling for 
the transverse gluons should slow down. It is even not inconceivable that 
the antiscreening running can be transformed into a screening one. In any 
case, one should expect that the value of the transverse coupling related to 
the matching scale rrid will be also essentially reduced in comparison with 
that in the vacuum QCD. Since the consideration in this section is rather 
qualitative, we adopt the simplest scenario with the value of the transverse 
coupling at the matching scale mj also coinciding with g s . 

The interaction potential between two static quarks in this theory at 
"short" distances r <~ m^ 1 reads: 

V(x, y, z) ~ 9 1 (14) 

where the dielectric constant e = 1 + ^ \ e q^\/ m q 1S ver Y large- Be- 

cause of the dielectric constant, this Coulomb like interaction is anisotropic 
in space: it is suppressed by a factor of y/e in the transverse directions 
compared to the interaction in the direction of the magnetic field. 

The potential (14) corresponds to the classical, tree, approximation 
which is good only in the region of distances much smaller than the con- 
finement radius tqc d ~ ^qcd ■ Deviations from this interaction are deter- 
mined by loop corrections. The analysis of the loop expansion leads to the 
following estimate of the new confinement scale Xqcd in QCD in a strong 
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magnetic field: 

/ \ b/b 

Xqcd " md [im ) ' (15) 

where b = (11N C - 2N f )/l2n and b = llN c /12n. Therefore, in a strong 
magnetic field, Xqcd is much less than Aqcd- 

The hierarchy Xqcd <C Aqcd is intimately connected with a somewhat 
puzzling point that the pairing dynamics decouples from the confinement 
dynamics despite it produces quark masses of order Aqcd or less [for a 
magnetic field all the way up to the order of (10 TeV) 2 ]. The point is that 
these masses are heavy in units of the new confinement scale Xqcd and the 
pairing dynamics is indeed weakly coupled. 



6. Arkady question: QCD with a large number of colors 

I did not discuss the case of QCD with a large number of colors in my talk. 
The question about it was raised by Arkady during the talk (he actually 
asked not myself, the speaker, but my coauthor Igor Shovkovy who was 
sitting next to him). Igor and myself got the answer the next day. 

Just a look at expression (7) for the gluon mass is enough to recognize 
that the dynamics in this limit is very different from that considered in the 
previous sections. Indeed, as is well known, the strong coupling constant a s 
is proportional to 1/N C in this limit. More precisely, it rescales as a s = ff-, 
where the new coupling constant a s remains finite as N c — > oo ('t Hooft 
limit). Then, expression (7) implies that the gluon mass goes to zero in this 
limit. This in turn implies that the appropriate approximation in this limit 
is not the improved rainbow approximation but the rainbow approximation 
itself, when both the vertex and the gluon propagator in the SD equation 
(4) are taken to be bare. This leads to the following expression for the 
dynamical mass of quarks [ 2] : 

1/2] 

/ " 1 v c \ 

-IT 



m 2 = C'\e q B\ cxp 



(iV2 - l)a t 



(16) 



where the constant C is of order one. The confinement scale Xqcd is close 
to Aqcd in this case [ 2]. 

It is natural to ask how large N c should be before the expression (16) 
becomes reliable. One can show [ 2] that the threshold value N* hr grows 
rapidly with the magnetic field [N l c hr > 100 for \eB\ > (lGcV) 2 ]. Expres- 
sion (16) for the quark mass is reliable for the values of N c of the order of 
jythr or \ aT g eT Decreasing N c below 7V* /lr , one comes to expression (1). 



12 



7. Conclusion 

QCD in a strong magnetic field yields an example of a rich, sophisticated 
and (that is very important) controllable dynamics. 
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